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Abstract. An abstract should be given 



abstract Almost 60 % of the B2 low luminosity ra- 
dio galaxies have been observed with the Hubble Space 
Telescope. We present an analysis of the dust features, 
which are often present in the form of circum-nuclear disks 
or lanes, and show that there are correlations between ra- 
dio source and dust properties. It is found that nearby ra- 
dio sources in which a jet has been detected tend to have 
dust more often than sources without jets; the dust is of- 
ten in the form of disks or lanes. Moreover the radio jets 

, are close to perpendicular to the disk or lane in the weaker 
radio sources (with P < 10^'' WHz~^). In stronger sources 
the orientation effect appears to be weak or even absent. 
Also the dust masses found in the weaker radio sources are 

, smaller than in the stronger ones (logM/MQ ~ 3 against 
5 respectively). More generally it appears that there is a 
correlation between dust mass and total radio power (for 
those sources in which dust has been detected); we show 
that this correlation is not induced by redshift. 

Key words. Galaxies: active; Galaxies: elliptical and 
lenticulars; Galaxies: nuclei 



1. Introduction 

A long standing question in radio astronomy has been why 
some elliptical galaxies host strong radio sources and oth- 
ers not. The reason may have something to do with the 
physical conditions near the active nucleus from which 
the radio source emanates. Although attempts have been 
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made in the past to link the presence of optical features 
like dust lanes and a radio source these have always been 
inconclusive due to the lack of adequate angular resolu- 
tion especially in the optical images. The quality of opti- 
cal imaging has now dramatically improved: imaging with 
the Hubble Space Telescope (HST) has provided valuable 
new information on the optical structure near the nucleus 
of radio galaxies. In particular, the HST/WFPC2 snap- 
shot survey of the SCR sample has produced a large and 
uniform database of images essential for a statistical anal- 
ysis of their host galaxies (Martel et al. 199£, De Koff et 
al. 1996). These studies revealed the presence of new and 



interesting features, some of them almost exclusively asso- 
ciated to low luminosity FR I radio galaxies. For example, 
the HST observations have shown the presence of dust 
in a large fraction of weak (FR I) radio galaxies which 
takes the form of ext ended nuclear disks (Jaffe et al 
De Koff et al. 



1993 



1996, De Juan, Colina & Golombek 



1996 



Verdoes Kleijn et al. 1999). Such structures have been 



naturally identified with the reservoir of material which 
will ultimately accrete into the central black hole. The 
symmetry axis of the nuclear disk may be a useful indica- 
tor of the the rotation axis of the central black hole (see 
e.g. Capetti & Celotti 1999| ), although it must be admit- 
ted that the precise relationship between these two axes 
remains uncertain. FR I sources represent an essential in- 
gredient in the Unified Schemes of radio-loud AGNs as the 
mis-oriented counter part of BL Lac objects (see Urry and 
Padovani 1995 for a review). Indeed, recent radio stud- 



ies have provided strong evidence in favour of relativistic 
beaming in their jets (Laing et al. 1999). Host galaxies 



of BL Lacs have been studied with the HST by Urry et 
al. ( 1999| ), and a similar analysis of galaxies hosting FR 
I sources should clarify if there is continuity in the prop- 
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erties of the respective nuclear regions, i.e. if any differ- 
ences are due solely to orientation effects. Furthermore, it 
was shown by Chiaberge, Capetti and Celotti (1999, 2000) 
that the optical nuclear sources commonly found in FR 1 
might represent the counterpart of the synchrotron radio 
cores. If this picture is correct, one expects that the dust 
features found in FR I sources may provide further indica- 
tions on the orientation of the radio source, provided that 
these features have something to do with the radio prop- 
erties. As many of the most extensively studied FR I radio 
galaxies are part of t he B2 sample of low luminosity radio 
galaxies (Fanti et al. 1987 ), this prompted us to perform 
a complete, high resolution optical study of these ^ 100 
sources. The radio characteristics of the B2 sample are 
different and in some sense complementary to those of the 
3CR already studied with HST and they can fill the gap 
between the " radio quiet & normal" and the " radio loud" 
ellipticals. Although only part (almost 60 %) of the B2 
sample has been observed with the HST, the high resolu- 
tion optical data can be considered as representative and 
unbiased with respect to the entire B2 sample. The new 
HST data were presented in Capetti et al. (2000). A statis- 
tical study of such a large sample of low luminosity radio 
galaxies enables us to establish how frequently we can de- 
tect optical jets, optical nuclear sources, circum-nuclear 
disks and dust lanes and if there is any relationship with 
the radio properties. The organization of this paper is as 
follows: in Sect. ||we briefly describe the HST data of the 
B2 sample and give some general information; in Sect. ^ we 
describe how we produced images in which the properties 
of the circum-nuclear dust can be analyzed; a discussion 
of the dust properties in given in Sect. ^. Throughout this 
paper we use a Hubble constant Hq = 100 km s^^Mpc^^ 
and qq ~ 0.5. 



2. HST observations of the the Sample 

A brief description of the B2 sample of low luminosity ra- 
dio galax ies and its present status was given in Capetti et 
al. ( 200C| ). The B2 sample of low luminosity radio galaxies 
has been extensively studied at radio wavelengths, espe- 
cial ly sin ce the 1980's (see Fant i et al. 19871 Parma et 
al. 



1987t De Ruiter et al. |l990t Morganti et al. |1997D 



Comparison with Einstein Satellite X-ray data was done 
by Morganti et al. (1988). A number of B2 radio galaxies 
were subsequently observed with ROSAT (Feretti et al. 



1995, Massaglia et al. 1996, Trussoni et al. 1997, Canosa 
et al. |1999[). 



Optical work on the entire B2 sample has somewhat 
lagged behind, but a complete broad band imaging survey 
of the B2 sample was carried out by Gonzalez-Serrano et 
al. ( |1993D and Gonzalez-Serrano & Carballo ( |200C| ), while 
narrow band Ha images were obtained by Morganti et al. 
( |l992|) . The IRAS prop erties of the sample were studied 
by Impey & Gregorini ( 1993 ). 

In the course of our HST program observations were 
obtained for 41 B2 galaxies, while public archive data exist 
for 16 additional objects, usually because these sources are 



also part of the 3C catalog (see Dc Koff ct al. 1996, Martel 
et al. |1999| , Verdoes Kleijn et al. |1999D . Therefore up till 
the present HST imaging has been done for 57 of the 100 
radio galaxies. In the following we will speak loosely of the 
B2 sample, when we discuss the sample of these 57 radio 
galaxies, while in the following the term 3C sample refers 
to 3C radio galaxies used by de Koff et al. (200C) that are 
not part of the B2 sample. 

Our program observations were obtained between 
April 9th and September 1st 1999, using the Wide Field 
and Planetary Camera 2 (WFPC2). The pixel size of 
the Planetary Camera, in which all targets are located, 
is 0.0455 arcsec and the 800 x 800 pixels cover a field 
of view of 36 x 36 arcsec. Two broad band filters were 
used, namely F555W and F814W, which cover the spec- 
tral regions 4500 - 6500 A and 7000 - 9500 A respectively. 
Although their transmission curves do not match exactly 
those of standard filters we will usually refer to them as 

V and / filters. The data have been processed through 
the PODPS (Post Observation Data Processing System) 
pipel ine for bias removal and fiat fielding (Biretta et al. 
|l996 ). One image was taken through each filter and the 
exposure time was always set to 300 s (with the excep- 
tion of B2 2116+26, which was observed for 350 s in the 

V band, as / band images were already available in the 
archive) . 



3. Circum-nuclear dust in low luminosity radio 
galaxies 

In order to study the dust in detail, we made use of the fact 
that we have both V and / images (at least for the sources 
that were observed in our program). Since the galaxy color 
V — I has only very small gradients the images in V and 
I are very similar and differ only by a constant. However, 
the dust features are more prominent in the V image, and 
this suggests we can fit the galaxy in the / image by a 
regular elliptical distribution and use this, to a first ap- 
proximation, as a dust-free model. Taking into account 
the constant offset between V and /, the V image was 
divided by the galaxy model as obtained from the / im- 
age. The resulting image of ratios will clearly show any 
dust features and thus can be easily analyzed; in regions 
without absorption the pixel values will be close to unity, 
while absorption results in values < 1. 

The dust features were classified as disks, lanes (pre- 
sumably edge-on disks) , filaments (irregularly curved fea- 
tures, possibly eccentric with respect to the galaxy nu- 
cleus), or complex (patchy structures). The classification 
of structures can be checked with the help of the album of 
images given in Capetti et al. ( ^000[ ). In Figs. |l| and ^ we 
give two of the best examples of disk-like structures found. 
These residual images (or "ratio maps") were analyzed us- 
ing the soft ware p ackage SYNAGE-I--I- recently developed 
by Murgia (2000). The arrows give the direction of the 
main jet. White colors represent absorption and therefore 
dust features. 
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Fig. 1. Residual image of B2 0034+25; the arrow indicates 
the direction of the main jet 
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Fig. 2. Residual image of B2 1322+36; the arrow indicates 
the direction of the main jet 

By no means all objects show absorption features: al- 
most half (42 %) of the observed galaxies have a very 
smooth light distribution without any sign of circum- 
nuclear dust. When dust is seen it is in most cases in the 
form of lanes or disks. Only in a minority of the B2 galax- 
ies do we see more complicated structures, like distorted 
filaments. A good example of c ompl ex dust structures is 
B2 1350+31 (see Capetti et al. |2000|) . 



4. Discussion 

4.1. The detection of dust 



We tried to make quantitative estimates of the amount of 
dust present, by considering a positive detection of dust 
only if in nine pixels that form a single area in the ra- 
tio maps the pixel value was less than 0.85, (see also de 
Koff et al. 2000). We calculated the masses of the clumpy 
dust features in the usual way, following Sadler & Gerhard 
(1985) and de Koff et al. (200C), by determining the cov- 
ering factor of the dust; this was done by summing all 
pixels in the ratio maps which had a value <0.85, such 
that the dust mass is: Afdust = S x A\ x F^^, where S is 
the area covered by the dust, Ax the mean absorption of 
the area under consideration and T the mass absorption 
coefficient, taken as 6 x 10^ mag kpc^Af-i (see de Koff 



et al. 2000). The mass is based on the absorption in the V 
images (see Sect. ^). The mean visual absorption is in the 
range 0.2-0.5 mag for all objects. The observed properties 
of the dust features (morphology, mean absorption and 
total area) are given in Table g. Sinc e 3C sources already 
had been analyzed by de Koff et al. ( 2000 ), we did not re- 
peat the calculation of A\, S and this is accounted for by 
an asterisk in the respective columns of Table |l|. A non- 
detection of dust means that less than nine pixels were 
found with value < 0.85 in a single area. Corresponding 
upper limits are based on this. In Table ^jwe list the up- 
per limits to the dust masses for those sources in which 
no dust was detected. 

We should remark that the detection of dust lanes or 
disks around galactic nuclei suffers from a number of se- 
lection effects, which are not so easily quantifiable. First, 
if the circum-nuclear dust takes, as we believe, the form 
of disks, then its detection will clearly depend on orien- 
tation of the disk with respect to the line of sight. Dust 
disks will be more easily seen if they are close to edge-on 
(in which case they will be called "lanes"); up till now 
the only counter example in the B2 sample is B2 0104+32 
(3C 31), in which we are able to see fairly complete el- 
lipses of dust around the nucleus. Another obvious selec- 
tion effect is due to redshift: since we require a minimum 
of nine pixels with significant absorption in a contiguous 
area, the detectable dust mass depends on redshift (see 
Sect. 4.3). This effect is shown in Fig. ^, where we plot 
the dust masses as derived for B2 sources, including the 
upper limits, as a function of redshift. Is should also be 
noted that the distribution of detected masses and upper 
limits appears to be discontinuous. 

Using the criterion for detection we find that dust is 
frequently present: 30/57 (53 %) of the galaxies show dust 
features, either in the form of bands or disk-like struc- 
tures, or more irregular patches. A summary of some gen- 
eral properties derived for these thirty sources is given in 
Table ||. Most columns of the table are obvious; it should 
only be mentioned that we give two determinations of dust 
masses; the first, in Col. 6, is the mass calculated from the 
absorption map as explained in this paper, while Col. 7 
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Table 2. B2 radio galaxies: sources without dust 



Name 


Redshift 


log P/WHz-^ 
1415 MHz 


log M/Mo) 
upper limit 


0055+26 


0.0472 


24.61 


2.3 


0120+33 


0.0164 


22.30 


1.5 


0708+32 


0.0672 


23.51 


2.6 


0722+30 


0.0191 


22.77 


1.6 


0755+37 


0.0413 


24.49 


2.2 


0924+30 


0.0266 


23.52 


1.9 


1003+26 


0.1165 


24.01 


3.0 


1005+28 


0.1476 


24.25 


3.2 


1101+38 


0.0300 


23.97 


2.0 


1113+24 


0.1021 


23.65 


2.9 


1204+34 


0.0788 


24.47 


2.7 


1251+27 


0.0857 


25.37 


2.8 


1257+28 


0.0239 


23.08 


1.8 


1422+26 


0.0370 


24.00 


2.1 


1450+28 


0.1265 


24.50 


3.1 


1502+26 


0.0540 


25.36 


2.4 




0.0931 


23.82 


2.8 


1521+28 


0.0825 


24.58 


2.8 


1553+24 


0.0426 


23.57 


2.2 


1557+26 


0.0442 


22.81 


2.3 


1610+29 


0.0313 


22.93 


2.0 


1613+27 


0.0647 


24.03 


2.6 


1615+32 


0.1520 


25.79 


3.2 


1658+30 


0.0351 


23.88 


2.1 


1827+32 


0.0659 


24.07 


2.6 


1833+32 


0.0586 


25.07 


2.5 


2236+35 


0.0277 


23.47 


1.9 



gives the dust mass (or an upper Umit) derived from the 
IRAS IR fluxes (if available). A "J" in Col. 8 means that 
the source has a radio jet. The percentage of dust detec- 
tions in the B2 sample may be slightly higher than in the 
3C sources, which in 33 % of the cases show dust (De Koff 
et al. 20001 ). However, it does not appear that there is any 
difference between FR I and FR II objects in the B2 sam- 
ple: the FR I sources with dust are 20/35 (57 %), against 
6/12 (50 %) of the FR I-II and FR II sources (for the FR 
classification see Capetti et al. 2000| ). Note that nine B2 
sources are either compact nuclei or non-classifiable. 

Our data strengthen the findings by Martel et al. 
( |1999| ) who noted that disk-like structures are preferen- 
tially associated to FR I type radio sources. Since FR I 
sources often have detected radio jets, it is expected that 
the presence of disklike structures is also related to the 
presence of jets. This is indeed the case, but we can even 
make a more detailed statement: if a source has complex 
dust features, i.e. no disk or lane, then most of the sources 
have no radio jet (7/9), while a majority of objects with 
a disk-like feature also have a radio jet (16/21). Another 
way of saying this is that the large majority (16/18) of the 
B2 galaxies with a radio jet also have a disk or lane (see 
Table Fisher's exact probability test (see e.g. Siegel 
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Fig. 3. Dust masses as a function of redshift; upper limits 
are based on the requirement that at least 9 pixels in a 
contiguous area should have a value <0.85 in the residual 
image 



1956) gives a probability of 0.9 % that the distribution in 
Table ^ is due to chance. 

It should also be mentioned that if we limit ourselves 
to sources with a jet and with redshift z < 0.03 (like the 
objects discussed by Verdoes Kleijn et al. 199E), then dust 
is detected even in 80 % of the cases, against the 53 % 
overall. This increase in the detection of dust features in 
more nearby galaxies regards only sources with a radio 
jet: if no jet is present the fraction of sources with dust is 
close to 50 %, both for nearby sources and for the sources 
with z > 0.03. 

We find no connection between dust content and radio 
size: the sources with dust have a median linear size of 
94 ± 35 kpc, against 76 ± 30 for the sources without dust. 

Table I gives the relation between the detection of dust 
and the presence of an optical core. The data on the opti- 
cal cores were taken from Capetti et al. (2002). It is quite 
obvious that there is no significant preference for sources 
with dust to have no detected optical core, or, reading 
the table differently, for galaxies without optical core to 
be more dusty. Although one would expect that the pos- 
sibility to detect an optical core is negatively influenced 
by the presence of dust, in line with the usual picture of 
the optical core being surrounded by a dusty disk or torus 
which obscures the core when the disk is seen close to 
edge-on, the absence of any significant effect may be in- 
fluenced by the strong selection effects discussed earlier: 
the detection of dust and the optical core in the standard 
model both depend on the orientation with respect to the 
observer. Especially the detection of dust may be sensi- 
tive, in a complicated way, to the redshift of the galaxy as 
well. However, perhaps most important of all, Chiaberge, 
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Table 1. B2 radio galaxies: properties of dust 



5 



Name 


Dust 


Major Axis 


Minor Axis 


r A (^QUSt j 


A\ 


V 

Zj 




type 


kpc 


kpc 


deg. 


mag 


pixels 


0034+25 


disk 


0.35 


0.12 


160 


0.37 


162 


0055+30 


disk 


0.24 


0.07 


43 


0.25 


210 


0104+32 


disk 


0.82 


0.65 


142 






0116+31 


filament 


5.02 


0.08 


20 


0.42 


2300 


0149+35 


disk 


0.84 


0.49 


170 


0.31 


5087 


0648+27 


disk 


0.28 


0.16 


45 


0.40 


1070 


0908+37 


disk 


0.31 


0.16 


135 


0.28 


31 


0915+32 


disk 


0.42 


0.34 


125 


0.37 


282 


1003+35 


disk 


2.85 


1.29 


60 






1217+29 


complex 


- 


_ 


- 


0.42 


3138 


1256+28 


disk 


0.10 


0.04 





0.27 


46 


1321+31 


disk 


0.23 


0.07 


43 


0.33 


710 


1322+36 


disk 


0.19 


0.05 


88 


0.39 


293 


1346+26 


filament 


4.02 


_ 


134 


0.35 


1307 


1339+26 


lane 


0.19 


0.05 


175 


0.28 


29 


1347+28 


lane 


0.36 


- 


140 


0.24 


9 


1350+31 


complex 


- 


_ 


- 






1357+28 


lane 


0.33 


0.08 


77 


0.23 


35 


1430+25 


disk 


0.21 


0.09 





0.28 


41 


1447+27 


filament 


1.18 


0.10 


138 


0.28 


604 


1455+28 


filament 


4.59 


1.68 


158 


0.48 


1881 


1457+29 


lane 


3.64 


0.80 


50 


0.43 


860 


1511+26 


lane 


0.30 




120 






1525+29 


disk 


0.43 




144 


0.24 


15 


1527+30 


disk 


0.48 


0.09 





0.34 


45 


1626+39 


filament 












1726+31 


lane 


9.00 




127 






2116+26 


disk 


0.44 


0.12 


65 


0.37 


10032 


2229+39 


disk 


1.30 


0.65 


169 






2335+26 


disk 


1.10 


0.76 


8 







Capetti & Celotti ( 1999 ) concluded that the optical nu- 
cleus is very often not obscured even if dust is present, so 
that the dust must be optically thin. This provides a nat- 
ural explanation for the absence of a correlation between 
dust and optical nucleus. In fact, all objects present a 
mean absorption consistent with an optical depth r < 1. 



4.2. Orientation of dust features and the direction of 
the radio axis 

Alignment effects involving the radio axis and some optical 
axis (for example the major axis of the galaxy) have been 
searched for in the past, but the results were mostly incon- 
clusive. However in nearby galaxies harbouring a (weak) 
radio source an alignment between the central dust fea- 
ture and the radio axis was found (Kotanyi & Ekers |1979D . 
Naively one would expect that the directions of the radio 
axis and the disk or lane are strongly correlated. Indeed, 
there is a general tendency for the radio axis to be perpen- 
dicular to the disk or lane, as can be seen in Fig. ^. There 
may also be a tendency for a more prominent alignment 



if a radio jet is detected (see Fig. ||), however, the effect, 
if it exists at all, is weak. 



According to de Koff et al. ( 1996 ), who used recent 
HST images of 3C sources, the alignment is most pro- 
nounced in sources with (small) disk structures. Figure || 
partly confirms this statement. Indeed the difference be- 
tween simple (i.e. disk-like) and more complex structures 
in non-B2 3C galaxies (lower panel) is significant at a 
confidence level of 98.5 %. However in our data on the 
B2 galaxies (upper panel) no significant difference is seen 
(confidence level < 90 %). The mean alignment angles 
are 62° for both simple and complex dust features in 
B2 sources, and 68° and 38° for simple and complex 3C 
sources respectively. In this and other figures we have de- 
fined disklike structures (including straight lanes) as sim- 
ple, and all other dust features (including distorted lanes) 
as complex. 

As we remarked in the previous Section, the simple 
structures tend to be found in sources with radio jets. It is 
well known that jets are more easily detected in weaker ra- 
dio sources of type FR I and it is therefore to be expected 
that dusty disks and lanes are more frequently present 
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Table 3. B2 radio galaxies: general properties 



Name 


z 


logP/W/fz"^ 
1415 MHz 


PA (radio) 
deg 


|APA| 
deg 


log M IMp. 
dust 


log M/Mp, 
IRAS 


Jet 


0034+25 


0.0321 


23.20 


93 


67 


3.6 


<6.3 


J 


0055+30 


0.0167 


24.08 


130 


87 


2.9 


5.9 


J 


0104+32 


0.0169 


24.21 


160 


18 


5.1 


6.4 


J 


0116+31 


0.0592 


24.95 


100 


80 


5.3 


7.0 


- 


0149+35 


0.0160 


22.33 


65 


75 


4.4 


6.2 


J 


0648+27 


0.0409 


23.62 


- 


- 


4.6 


7.3 


- 


0908+37 


0.1040 


24.84 


15 


60 


3.7 


<7.2 


J 


0915+32 


0.0620 


24.00 


30 


85 


4.4 


<6.8 


J 


1003+35 


0.0989 


25.78 


120 


60 


6.3 


<6.6 


J 


1217+29 


0.0021 


21.24 


- 


- 


2.6 


4.7 


J 


1256+28 


0.0224 


23.05 


- 


- 


2.6 


<5.4 


J 


1321+31 


0.0161 


23.85 


120 


77 


3.6 


<5.7 


J 


1322+36 


0.0175 


24.55 


15 


73 


3.3 


<5.8 


J 


1346+26 


0.0633 


24.55 


25 


71 


5.0 


<6.5 


- 


1339+26 


0.0757 


24.30 


30 


35 


3.4 


<7.0 


- 


1347+28 


0.0724 


24.05 


50 


90 


2.7 


<7.0 


J 


1350+31 


0.0452 


25.03 


90 


- 


6.1 


6.8 


J 


1357+28 


0.0629 


24.03 





77 


4.4 


<6.8 


J 


1430+25 


0.0813 


24.20 


20 


20 


3.6 




- 


1447-1-27 


fl 0306 


22.78 


27 


69 


4.1 






1455+28 


0.1411 


25.22 


38 


60 


5.9 






1457+29 


0.1470 


24.89 


155 


75 


5.5 






1511+26 


0.1078 


25.34 


20 


80 


4.9 






1525+29 


0.0653 


23.98 


30 


66 


3.0 


<6.8 


J 


1527+30 


0.1143 


24.05 


45 


45 


4.0 


<7.3 




1626+39 


0.0303 


24.49 


90 




4.9 


<5.7 


J 


1726+31 


0.1670 


25.89 


110 


17 


5.6 


<7.6 




2116+26 


0.0164 


22.79 


157 


88 


4.7 


6.3 


J 


2229+39 


0.0181 


24.03 


10 


21 


4.3 


6.3 


J 


2335+26 


0.0301 


24.88 


125 


63 


3.8 


<6.0 


J 



Table 4. Detection of dust and radio jets 





Radio sources 
with jet 


Radio sources 
without jet 


Total 


Galaxies with dust-disk or lane 


16 


5 


21 


Galaxies with complex dust 


2 


7 


9 


Total 


18 


12 


30 



in the weaker radio sources with jets; any alignment be- 
tween radio axis and dust-disks should be sought there. 
This is illustrated in Fig. ^ where we plot the alignment 
angle against radio power. Clearly many sources, and ac- 
tually all of the weak ones with P < 10^'*WHz~^, have an 
alignment angle > 60°. For the FR II sources (coinciding 
mostly with the non-B2 3C sources in Fig. ^ there are 
quite a few sources with small angles (< 60°). The mean 
alignment angle is 75.3 ± 8.7° for the B2 sources with 
P < lO^^WHz-i, against 57.3 ± 25° for the stronger B2 
sources. The difference between the average angles is not 
significant, but the variance of the distributions is differ- 



ent (confidence level > 99 % from Fisher's F-distribution). 
This of course reflects the fact that all the weak sources 
are clustered towards high alignment angles. 

4.3. Dust masses 

As explained in Sect. 4.1, we calculated the masses of the 
dust features in the usual way, following Sadler & Gerhard 
(1985). In Fig. ||we show the derived dust masses against 



the alignment angle. No correlation is seen. 

A confirmation that the simple structures have on av- 
erage lower masses is shown in Fig. ^: the lower masses 
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Table 5. Detection of dust and optical core 



7 





Galaxies with 


Galaxies without 


Total 




optical core 


optical core 




Galaxies with dust 


9 


21 


30 


Galaxies without dust 


9 


18 


27 


Total 


18 


39 


57 



12 



10 
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Fig. 4. Radio-dust alignments in B2 sources (horizontally 
shaded) and 3C sources that are not part of the B2 sample 
(diagonally shaded) 
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Fig. 6. Alignments according to dust structure. For defini- 
tion of simple and complex dust structure see text. Upper 
left: simple B2 sources; upper right: complex B2; lower 
left: simple 3C; lower right: complex 3C 
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Fig. 5. Radio-dust alignments according to the presence 
(upper panel) or absence (lower panel) of a radio jet 
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Fig. 7. Alignments according to radio power. Filled sym- 
bol: B2; open symbol: 3C sources 
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correspond to the disk and lane structures. For the B2 
galaxies with simple dust structures we find an aver- 
age < log(M/M0 >= 3.9, with a standard deviation of 
the mean of 0.20, while complex structures have an av- 
erage of 5.0 with standard deviation 0.36. According to 
Student's t-test this difference is significant at the 99.5 % 
level. A similar difference is obtained for the 3C galax- 
ies: < log(M/M0 >= 4.68 ± 0.25 for simple dust and 
5.2 ± 0.13 for complex dust. This difference is significant 
at a confidence level of 98 %. The higher mass values of 
the 3C galaxies reflect the fact that they are generally 
stronger radio sources. Anyway, it is possible that the ob- 
jects classified as containing complex dust structures also 
have small disks, hidden in the complexity of the dust 
features. 

Most striking, however, is the correlation shown in Fig. 
[lO| , between the total radio power and the detected dust 
masses. A linear regression between log P and log M/Mq 
gives a slope of 0.61 ±0.16. We can re ject (using the recipe 
given in e.g. Haber & Runyon 1973) the hypothesis that 
radio power and dust mass are uncorrelated at the 99.9 
% level. However, it is known that radio power and red- 
shift are correlated in flux limited samples; using the test 
discussed by Macklin ( |1982| ), we can reject the hypothe- 
sis that the apparent correlation between power and dust 
mass is due exclusively to a correlation between P and 
z and M and z, at the level of 4.3cr. We may illustrate 
this strong correlation (independent of redshift), by di- 
viding our sample in four redshift intervals (z < 0.02, 
0.02 < z < 0.05, 0.05 < z < 0.08 and z > 0.08). The 
correlation remains in three of the four redshift interval 
(see Fig. |l^); although it is formally not significant in the 
lowest redshift bin (0.45(7) the correlation is very strong 
indeed in the other redshift bins, even if the number of 
objects involved is only of the order of 7 or 8 per bin: the 
significance levels are 3.1, 4.3 and 4.9 fj in order of increas- 
ing redshift. The location of the upper limits in Fig. ^ 
strongly suggests that the distribution in the log P — log AI 
plane is bimodal, that is, if there is dust detected in the 
form of disks or lanes, then the amount of dust depends on 
radio power (or vice- versa), but in addition to this there 
are radio sources in which no dust is detected, practically 
independently of the radio power (the slow rising of the 
upper limits in Fig. |l^ with radio power is due to increas- 
ing redshift). 

No correlation is found between detected dust masses 
and the absolute magnitude of the galaxy and this is in 
line with the result of de Koff et al. ( ^000 ); therefore the 
relation dust mass against power stands on its own, and 
does not appear to be induced by other possible correla- 
tions between galaxy properties. 

In order to compare the correlation derived from the 
B2 sample with other data we show in Fig. 
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the same 

plot as in Fig. |10|, but now with non-B2 3C sources added 
and also an optically selected sample of nearby (z < .011) 
E and SO galaxies described in Tran et al. ( ^001 ). For 
clarity no upper limits were plotted. The galaxies in which 
both radio emission (based on the NVSS) and dust were 
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Fig. 8. Alignments according to dustmass. Filled symbol: 
B2; open symbol: 3C sources 
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Fig. 9. Masses according to dust structure (see text). 
Upper panels: B2, lower panels: 3C sources. Left: simple 
structures; right: complex structures 



detect ed are 18 of the 107 galaxies discussed by Tra n et al . 
( 2001 ). Since the selection criteria of the Tran et al. ( 2001 ) 
are rather complicated no statistical conclusions should 
be drawn from the location of these galaxies in Fig. ^ 
we just remark that they roughly fall in the region were 
we would expect them to be according to the correlation 
based on the B2 sample. 

In a number of cases we can also derive dust masses 
from IRAS data, although most B2 sources have only up- 
per limits. We u se the 100 ^m data given by Impey & 
Gregorini (1993), if available. The dust masses were cal- 
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Fig. 10. Dust mass as a function of total radio power (B2 
sources). Different symbols indicate different redshift in- 
tervals. Open circles: z < 0.02; open stars: 0.02 < z < 
0.05; filled squares: 0.05 < z < 0.08; crosses: z > 0.08 
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Fig. 11. Dust masses as a function of total radio power, 
for B2 (filled dots), 3C (open circles) and sources from 
Tran et al. (stars) 



culated using a dust temperature Td = 30 K (see de Koff 
et al. 200C| ); the results are shown in Fig. |l2|. The derived 
IRAS masses are typically a factor '^100 higher than the 
masses derived from our absorption analysis. To bring the 
masses down to the level of the masses determined here 
from the absorption analysis, we would need to increase 
the dust temperature in most cases to unrealisticly high 
values (» 100 K). Only for B2 0104+32 (3C 31), B2 
1350+31 and B2 2116+26 can we bring the dust masses 




log (Mr^Mj 

Fig. 12. The relation between dust masses derived from 
absorption, as calculated here and from IRAS observations 



in agreement with Td < 100 K. We therefore conclude that 
the detected IRAS dust masses are in general one or two 
orders of magnitude higher than the masses derived from 
the absorption analysis. The upper limits in Fig. ^ are 
consistent with this. It is of course not surprising that the 
IRAS masses are larger, as the IRAS fluxes measure the 
overall presence of dust (including diffuse components), 
while the absorption analysis is only sensitive to sharp and 
well defined features, which in many cases means only the 
dusty disk close to the nucleus. Apparently the circum- 
nuclear features measured by us represent not much more 
than about between 1 and 10 % of the overall dust. 



5. Summary 

1. Dust features are detected in slightly over 50 % (30/57) 
of the B2 sources, which is a bit higher than was found 
for 3C sources by de Koff et al. (|200(]| ). Ahhough this 
would suggest a difference between FRI (majority of 
B2) and FRII (majority in 3C) sources no significant 
effect is found. 

2. If a B2 radio source has jets, then in 16/18 sources the 
dust around the galactic nucleus comes in the form of 
disks or lanes, i.e. of simples structures; this correlation 
is very strong: the probability that the distribution is 
only due to chance is 0.9 % . Conversely, in seven of 
the nine sources with complex dust structures no radio 
jet is seen. 

3. The radio jets tend to be perpendicular to the dust 
lanes or disks, in particular the weaker FR I radio 
sources, with P < lO^^'WHz^^. Moreover small dust 
masses (with log M/Mq < 3) are found mostly in 
these weaker sources. In part this may be a selection 
effect, because the weakest sources are also the most 
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nearby, in which small dust masses are more easily de- 
tected. 

4. We have shown that if circum- nuclear dust is present 
in low luminosity radio sources, then its amount and 
its morphology are linked to the properties of the radio 
sources. In other words the radio source "knows" about 
the presence of dust. We base this conclusion on the 
fact that there is a well established correlation between 
the dust mass and the total (or core) radio power, as 
described in the previous Section. Since the radio and 
optical properties are (in observational sense) only in- 
directly dependent via the distance, the only way to 
explain this correlation away is through a selection ef- 
fect that works both on the derived dust masses and on 
the radio power. The fact that the correlation between 
radio power and dust mass remains valid also within 
different redshift intervals strongly suggests that the 
correlation is intrinsic. However this correlation holds 
only if a significant amount of dust is present, and there 
may be a discontinuity (or bimodality in the mass dis- 
tribution), as is also clear from Fig. |l0|. 
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